Witnessing Quantum Coherence: from solid-state to biological systems 
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Quantum coherence is one of the primary non-classical features of quantum systems. While 
protocols such as the Leggett-Garg inequality (LGI) and quantum tomography can be used to test 
for the existence of quantum coherence and dynamics in a given system, unambiguously detecting 
inherent "quantumness" still faces serious obstacles in terms of experimental feasibility and efficiency, 
particularly in complex systems. Here we introduce two "quantum witnesses" to efficiently verify 
quantum coherence and dynamics in the time domain, without the expense and burden of non- 
invasive measurements or full tomographic processes. Using several physical examples, including 
quantum transport in solid-state nanostructures and in biological organisms, we show that these 
quantum witnesses are robust and have a much finer resolution in their detection window than the 
LGI has. These robust quantum indicators may assist in reducing the experimental overhead in 
unambiguously verifying quantum coherence in complex systems. 



PACS numbers: 

Introduction 

Quantum coherence, or superposition, between differ- 
ent states is one of the main features of quantum systems. 
This distinctive property, coherence, ultimately leads to 
a variety of other phenomena, e.g., entanglement [TJ |2]. 
It is also thought to be the power behind several "quan- 
tum tools" , including quantum information processing[3], 
metrology [3] , transport |5j , and recently, some functions 
in biological organisms [B] (e.g., efficient energy trans- 
port). 

Identifying quantum coherence and dynamics in an ef- 
ficient way, given limited system access, is indispensable 
for ensuring reliable quantum applications in a variety of 
contexts. Furthermore, the question of whether quantum 
coherence can really exist in biological organisms in vivo, 
e.g., in a photosynthetic complex or in an avian chemi- 
cal compass, surrounded by a hot and wet environment, 
has triggered a surge of interest into the relationship be- 
tween quantum coherence and biological function [8] . 
In these cases, full-system access is often very limited, 
and signatures of quantum coherence are often indirect. 

The existing methods for identifying quantum coherent 
behavior can be generally classified into two types. The 
first type are based on imposing what can be thought 
of as a classical constraint [9], such as macroscopic re- 
alism and non-invasive measurements in the Leggett- 
Garg inequality (LGI) [ID] , or realism and locality in 
Bell's Inequality. Even though inequalities like the LGI 
were originally envisaged as a fundamental test of phys- 
ical theories, a violation of the LGI can also be con- 
sidered as a tool for classifying the behavior observed 
in experiment as quantum or classical. However, the 



Leggett-Garg inequality faces severe experimental diffi- 
culties when used as such a tool as it requires nonin- 
vasive measurements, e.g., via quantum nondemolition 
(QND) measurement [HHH], weak measurement [13], or 
quantum-gate-assisted ideal non-invasive measurements 
[14] , Because of this only a few tests of the LGI have 
been reported [T4TH5] . 

The second type of test is based on deduction; do the 
results of a given experiment sufficiently correspond to 
the predictions of quantum theory (or classical theory, 
depending on the approach). Quantum witnesses can be 
considered as one such test, as they use the knowledge 
of a quantum state or of some quantum dynamics to de- 
termine whether an experimental system possesses quan- 
tum properties. Some examples that have been employed 
elsewhere include witnesses of entanglement [THl dDJ > di- 
rect measurement of coherence terms of density matri- 
ces, or the analysis of process tomography [21] for non- 
classical state evolutions. The experimental realization 
of this kind of verification usually needs tomographic 
techniques, and then the required experimental resources 
in terms of measurement settings increases exponentially 
with the system complexity [T9] [21] [22] . Moreover, quan- 
tum state and process tomography are still difficult to 
implement in general systems and for general state evo- 
lutions, e.g., particularly in systems like charge transport 
through nanostructures, the transfer of electronic excita- 
tions in a photosynthetic complex, or systems where the 
state space is large. 

In this work, we introduce two quantum witnesses to 
verify quantum coherence and dynamics in the time do- 
main, both of which have various advantages and disad- 
vantages. Both are efficient in the sense that there is 
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no need to perform noninvasive measurements or to use 
quantum tomography, dramatically reducing the over- 
head and complexity of unambiguous experimental veri- 
fication of quantum phenomena. 

We apply these quantum witnesses to five exam- 
ples: (1) electron-pair tunnelling in a Cooper-pair box 
and coherent evolution of single-transmon qubit, (2) 
charge transport through double quantum dots, (3) 
non-equilibrium energy transfer in the photosynthctic 
pigment-protein complex, (4) vacuum Rabi oscillation 
in lossy cavities, and (5) coherent rotations of photonic 
qubits. Furthermore, as we will illustrate in these ex- 
amples, our quantum witnesses possess a finer detection 
resolution than the LGI. 

Both witnesses, which we will introduce shortly, in- 
volve the following steps: (Figure la): first, we prepare 
the system in a known product state with its environ- 
ment (or reservoir, here we use both terms interchange- 
ably) psr{0)- We then let psb,(0) evolve for a period of 
time to j to reach the state psii(to) (during which time one 
hopes the state has acquired significant coherence due to 
its internal dynamics). The second step is to implement a 
quantum witness using a "correlation check" between the 
state psii(to) and its state at another later time t > t , 
Psii(t)- The goal of this correlation check is to investi- 
gate non-classical properties in these two-time state-state 
correlations (see Figure lb). If the state psR(io) can be 
detected then by our quantum witness as having quan- 
tum properties, this implies that either the system state 
Ps{to) = Tr^[p5fl(to)] possesses significant quantum co- 
herence or that the state psn(to) is an entangled system- 
bath state. 



Results 

In order to find a signature of quantum dynamics we 
start by seeking characteristic features of classical dy- 
namics or states [23] • All separable mixtures of system- 
reservoir states, with no coherent components, which we 
call classical states, obey the following relation for their 
two-time correlations: 



a State Preparation b Correlation Check 



(Qm(t)Qn{to))n = Pn(to)^mn(t, t ). 



(1) 



See Methods for the proof. Succinctly put, equation ([!]) 
implies it is possible to define all future behavior based 
on only the system's instantaneous expectation values 
p n (to). However, most quantum correlation functions 
also obey this relation under certain measurement con- 
ditions. For example, a correlation function constructed 
from two-time projective measurements has this form as 
the measurement at to destroys the coherence in the state 
at that time. Here Qi is an observable which measures 
if the system is in the state i. This state is assumed 
to have a classical meaning (e.g., localized charge state, 
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FIG. 1: Detecting quantum coherence and dynamics. 

Generic procedure for detecting quantum coherence: We need 

(a) state preparation and (b) a correlation check. With the 
freedom to manipulate the system state, the initial state of the 
total system can be reasonably prepared as a product state 
Psr(0) — ps(0) CB> pi?(0), where pn(0) is the reservoir state. 

(b) and (c) show the correlation check, or measurements, we 
base our quantum witnesses on. Assuming the system state 
at time to is ps(to) = Qn(to = Pn{to)), the probability of be- 
ing in a state Q m at later time t is determined by the n — > rn 
propagator Q mn (t, to)- The general quantum correlator is de- 
fined by {Qm(t)Q n (to)) q (upper green path) and the classi- 
cal one is defined by p n {to)Qmn(t, to) (lower brown track). In 
this expression, Q mn (t, to) is the probability of measuring the 
system in state m at time t given that it was in the state 
n at time to- Both definitions describe the connections be- 
tween Q n (to) and Q m (t) for arbitrary states ps{to) with a 
distribution of state populations {p n (to)}- As shown in the 
Methods, all classical dynamics should satisfy the relation, 
{Qm(t)Q n (to)) Q = Pn(io)fimn(Mo), whereas violations of the 
equality reveal a signature of quantum dynamics. While our 
quantum witnesses are derived from this "correlation check" , 
the experimental requirements for each witness differ. 



etc) and the observable is normalized so that its expecta- 
tion value is directly equal to the probability of observ- 
ing the system in that state (Qi) = pi- The propaga- 
tor Q mn (t,to) is the probability of measuring the system 
in state m at time t given that it was in the state n 
at time t (and which in principle depends on the state 
of the reservoir, so can include classical non-Markovian 
correlations, see Methods). Several other recent tests of 
quantumness [HI [24 27 rely on imposing Markovianity 
on fl nm (t,to). In our first witness we avoid taking that 
approach so that we can still distinguish quantum from 
classical non-Markovian dynamics. However we will use 
it in our second witness. 

In principle, one could use Eq. Q to construct a quan- 
tum witness of the form: 



(Qm{t)Qn(to))n ~ Pn(to)Qmn{t, to) 



(2) 
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Where a non-zero result Wqq > 0, implies the state at to 
can be considered as quantum in that it contains quan- 
tum coherence which effect its future evolution. However, 
as mentioned above, most quantum correlation functions 
also obey equation which will give Wqq = 0. Is it 
ever possible to observe a non-zero Wqq? In some cases 
coherence, or "amplitude" , sensitive correlation functions 
are encountered in quantum optics [28], and in linear- 
response theory [53]. However, these are typically ex- 
tracted from spectral functions in the steady state, or 
put in a symmetrized form, in which case any affect on 
the correlation function from the initial state coherence 
may be lost. In all the examples we consider in this work 
this witness Wqq cannot be directly measured, as the 
initial coherence is of course destroyed by the first (pro- 
jective) measurement. Fortunately, Wqq, via Eq. 0, 
gives us a way to develop a more generally applicable 
and valid witness. 



Witness 1 

Our first practical witness (which is the main result 
of this work) can be derived from Eq. ([I]) by including 
normalization. Noting that all classical system-reservoir 
states obey, 



(Qm(t)) = J2pn(to)Cl mn (t,t ), (3) 
n=l 

where d is the number of states n in, or dimensionality 
of, the system state space, we define our first quantum 
witness as 



Wq := 



(Qm(t)) -53l»n(*o) n 

r; 



* (Mo! 



(4) 



If Wq > 0, we can define the state at to as quantum. 
Compared with the witness Wqq and the tests of the 
LGI, Wq can always be directly measured, and ideal non- 
invasive measurements are not necessary. In experimen- 
tal realizations, measuring the population-related quan- 
tities, or expectation values, (Q m (t)) and {p n {to)}, is 
generally more feasible than constructing full correlation 
functions, particularly in systems which rely on destruc- 
tive (e.g., fluorescence) measurements. Where correla- 
tion functions can be measured with projective measure- 
ments, the second term can of course be replaced with 

2 n Pn(*o)^ron(*,*o) = J2n (Qm (t)Qn (*o)) ■ 

However, determining all the propagators Q mn (t,to) 
with which to construct the witness requires, in principle, 
that we can prepare the system in each one of it states 
n exactly (or, alternatively if correlation functions con- 
structed from projective measurements are available, it 
requires that we measure every possible cross-correlation 
J2n(Qm{t)Qn(to))) ■ In the former case (where we use 



state preparation) we trade-off the need to do non- 
invasive state measurement with the need to perform 
ideal state preparation. In complex systems it may be 
difficult to prepare the system in each one of its states to 
construct these propagators, and in some cases we may 
not even have knowledge of the full state-space of the 
system. 

Importantly, this problem can be easily overcome by 
noticing that the individual terms in the sum in Eq. @ 
are always positive. Thus when constructing the sum we 
can stop as soon as the witness is violated by this par- 
tial summation (i.e., when the terms in the summation 
together are larger than (Q m (t))), reducing the experi- 
mental overhead substantially (see Figure 4 for a practi- 
cal example, where we show it is sufficient to include just 
one term in the sum of Eq. Q). 

Note that with this witness we do not distinguish 
between just system-coherence or quantum correlations 
(entanglement) between system and bath/reservoir (see 
Methods). In addition, if there are classical correla- 
tions between system and reservoir, i.e., classical non- 
Markovian effects [3D], then some additional experimen- 
tal overhead is needed to eliminate this from giving a 
"false positive" . If this overhead is ignored this represents 
a "loop-hole" in this witness, and in some situations may 
be an obstacle for its unambiguous application. We will 
discuss this explicitly later with an example of a pho- 
tosynthctic light-harvesting complex where the system 
and reservoir are strongly correlated both classically and 
quantum mechanically. 



Witness 2 

For our second witness we impose the extra condition 
that fi mn (t,*o) = Qmn{t' ,t'o) for t - to = t? - t' = T, 
for any time interval r. This assumption restricts us to 
a widely-studied subset of quantum processes where the 
system-bath/reservoir interaction is Markovian. We will 
show that, under the assumption that our system lies 
within this subset, quantum properties can be identified 
without needing to explicitly measure propagators (i.e., 
neither exact state initialization or non-invasive measure- 
ments are required). The trade-off in this case is that 
the witness cannot distinguish certain types of classi- 
cal dynamics (e.g., classical non-Markovian) , from quan- 
tum properties of the system. Still, this witness exceeds 
the tests proposed in earlier works under the same con- 
straints which still required either non-invasive measure- 
ments or state preparation [IS] [35] . 

This subset of quantum processes can be described as 
having weak coupling between system and reservoir so 
that system-reservoir state is always a product state, 
and the bath/reservoir state does not evolve in time, 
i.e., pn(t) — pr(0). A large number of systems exist 
in this regime [30], with well-developed models such as 



4 



the master equation under the Born approximation oper- 
ating within this class (see, e.g., [30-32] ) . For such cases, 
we can extend the first witness so that we replace the 
need to prepare the system state with that of needing to 
repeatedly measure expectation values (not correlation 
functions) a number of times that scales linearly with 
system size. To show this, we consider an extension of 
Eq. ([3]) involving a system of d linear equations repre- 
sented in matrix multiplication form as follows: 



where the d x d matrix Pj 
and £l rn j and Q T 



mj 



has elements [P 7 *]fe f 



(5) 



vectors with elements 



are d x 1 column 
ft m j]ni = n m „y](r) and 
(Qm(t\j,k]))> respectively. Here, t y jfc] 
and tu fci constitute the jth nontrivial time-domain set 



[Q 



mj \ k 1 



Tj-{to\j,k]tt\j,k]\ t\j,k] -to{j,k] = T ',k — 1, 2, ...,ri}. For a 
given time difference r and a time pair (toy,*] > t{j,k] ) G Tj, 
one can use the most experimentally-feasible method of 
measurement, i.e., invasive measurement, to obtain the 
information about the state populations, Pn(^o[j',fel) an d 
the expectation values (Q m (% w)). 
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FIG. 2: Detecting quantum oscillations, (a) shows a 
schematic circuit of a single-Cooper-pair box [361 137] . Its 
Hamiltonian is described by Hcb = Ec{l — 2n 9 )/2(|2)(2 — 
1 1)(1| ) - Ej/2(\ 1)(2| + 12)(1 1), where E j and E c are the Joseph- 
son energy and the single-Cooper-pair charging energy of the 
box, respectively. The relative energy of the state with no 
excess Cooper pairs in the box |1) and the state with one ex- 
cess Cooper pair |2), is controlled through the gate voltage, 
which is parametrized by n g . The resonance of the states 
|1) and 1 2} can be brought from the initial state |1) by the 
applied voltage pulse for n g — 0.5. (b), Detecting quantum 
dynamics in the resonance of the two charge states (n g = 0.5) 
with the quantum witness Wrm ■ We use the realistic param- 
eter Ej = 51.8 jj,eV. The positive regions are identified as the 
quantum areas. 

Given a set of measurement results to sufficiently de- 
scribe the state populations, the vector fl m j can be deter- 
mined by simple algebraic methods. For nonzero deter- 



minant det(Pj), we have ft mn ui(r) 



where P~] is the matrix formed by replacing the nth col- 



det(Pg)/det(P,), 



umn of P j by Q m j . For an arbitrary pair of time-domain 
sets, say Tj and Tji, we impose an additional condition 
(not used in the earlier witnesses) that their propagators 
should be identical for all classical systems (within the 
subset described above): fi mj - = O mj -/. If the system 
and its environment are classically-correlated, i.e., they 
are not in a product state, this assumption does not hold. 
Any comparison between Q m j and £l m j> can be consid- 
ered as a quantum witness for this subset, such as the 
vector-element comparison: 



det(Pg)det(P J .) - det(Pg,)det(P j ) 



(6) 



raj 



If yVsimn > 0, and under the assumptions described ear- 
lier, we can again assume that some of the (set) of initial 
states are quantum. Since measuring Wsimn requires the 
information about state populations only and can be per- 
formed with invasive observations, implementing Wnmn 
can be more practical than implementing Wqq (|2| and 



Examples 

To illustrate the effectiveness of our witnesses we now 
present five example systems where they could be ap- 
plied. For each example we choose which ever witness is 
more appropriate, given the properties of that system. 

Rabi oscillations in superconducting qubits 

The oscillations of state populations are commonly 
thought of as a signature of quantum dynamics. The 
measurement of these kind of oscillations is widely em- 
ployed for many experiments. The observation of such os- 
cillations alone, however, is not definitive evidence for the 
existence of quantum coherent dynamics and can even be 
mimicked by the solutions of classical autonomous rate 
equations, e.g., Ref. [531 154] , 

As a first example of the application of our witnesses 
we apply Wa mn to a two-level system composed 
of the two lowest-energy states in a single-Cooper-pair 
box [3"5rf3"T] . Figure 2a. We can take n = 1, m = 2, 
for example together with the designation Tf. {trjufc] = 
(*+J-l)*o,t[i,fc] = (k+j-l)t +T\k = 1,2} for j = 1,2, 
Figure 2b illustrates that the quantum witness WV221 de- 
tects the presence of quantumness in the Cooper-pair 
tunneling. Since only information about state popula- 
tions is required, this witness is easy to apply in practice 
with simple invasive measurements and can be readily 
applied to the existing experiments in the time domain 
|36| 137] without any additional experimental overhead. 

One can also consider an application of our witnesses 
to single- and multiple-transmon qubits coupled to trans- 
mission lines in circuit quantum electrodynamics [381139] . 
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FIG. 4: Detecting quantum properties of the FMO 
complex. Magnitude of the first witness Wq [Eq. Q for 
t — to = t] for the FMO complex assuming the final measure- 
ment is done on site m = 1 of the seven site FMO complex, 
for both (a) T = 77 K and (b) T = 300 K. A detection 
via our first witness is clearly visible for an initial evolution 
greater than to = 0.3 ps at 77 K. In comparison, for same pa- 
rameters we employ for the witness, the LG inequality only 
reveals a violation for upto 0.035 ps 44]. (c) and (d) show 
the first witness with only limited access, i.e. with only state 
preparation and measurement on site 1. Quantum coherence 
is only detected when pi(to)r2n(t, to) > (Qi(*)}- In all figures 
the bath parameters used were 7 -1 = 50 fs and A = 35 cm -1 , 
and the Hamiltonian is the same as that used in Ref . [8] . For 
the Hierarchy calculation, we used the "Ishizaki-Tanimura" 
truncation scheme and truncation as taken at K = and 
N c = 8 (see Methods, or Ref. [5], for the meaning of these 
parameters). 



FIG. 3: Detecting quantum transport through a dou- 
ble quantum dot. a, Schematic of a single-electron dou- 
ble quantum dot (DQD). Here we assume that the DQD is 
weakly coupled to leads under a large bias. Its Hamilto- 
nian is Jfcqo = A(\L){R\ + \R){L\) with the electron state 
basis {|L),|i?},|0}} where A is the tunnelling amplitude be- 
tween the left-dot and right-dot electron states \L) ,\R). The 
transport between dots and leads is described by the self- 
energy, E[p] = -l/2E a=t ,ji r »[ s « s ip ~ 2sips a + ps a si], 
where sl = \0){L\, sr = |-R}(0|, and Fl and F_r are the left 
and right tunnelling rates, respectively. We assume charge 
detectors (CDs) are used for the measurements, but invasive 
current measurements are also sufficient (not shown here). 
b,c, Verifying quantum transport through DQD with Wq 
[Eq. |E| for t - 1 — r] and Wn 32 [Eq. , respectively. Here 
we define |0), \L), and \R) by |1), |2), and |3), respectively. 
For the setting Tl = 4, Yr — 0.1, and A = 1, the non- 
vanished Wq and cWn32 indicate the quantum-transport re- 
gions, where c = (piP2f>3) _1 for the stationary state. 



where qubit-state measurements are performed by mon- 
itoring the transmission through the microwave cavity 
[39] ■ For the simplest case of one-qubit rotation, the co- 
herent evolution is driven by the Hamiltonian[38 



ff = ^|l)(l|+ e (t)(|0)(l| + |l)(0|), 



(7) 



where e(t) is the microwave pulse to induce transitions 
between qubit states |0) and |1) with an energy differ- 
ence Hui. Through properly choosing the pulse e(t), a 
reliable single-qubit gate, e.g., the Hadamard transfor- 
mation (H), can be created. Here, we use the quantum- 
process-tomography-based optimal control theory [4(J to 
design the microwave pulse for such a gate (£n) with a 
process fidelity of about 94%. We use the first witness 
Wq in the form: 



(0(£ 2H )> - 5>n(£ H )Oo n (£ H ; 



(8) 
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to show that the process £h creates coherent rotations. 
When setting the input state as |0), the value of our wit- 
ness is about Wq « 0.45, which certifies the quantumness 
of £ H - 

Quantum transport in quantum dots 

Experimentally distinguishing quantum from classical 
transport through nanostructure remains a critical chal- 
lenge in studying transport phenomena and designing 
quantum electronic devices. As mentioned in the in- 
troduction, using time-domain methods to verify quan- 
tum coherence, such as by testing the Leggett-Garg in- 
equality, can be very demanding. We illustrate here how 
our witnesses are valid under invasive measurements by 
modelling single-electron transport through double quan- 
tum dots (Figure 3a). Compared with the time periods 
identified by the Leggett-Garg-type approach [35], the 
quantum witnesses Wq (Figure 3b) and Wn mn (Figure 
3c) can detect a much larger quantum coherence win- 
dow. For Wnmn, we employ the settings Tf. {tow.fcl = 
[k + d(3-l)]t ,t {l ,k] = [k + c'(i - l)]to + r\k = 1,2,3} 
for j = 1, 2. Here d is large such that the whole system 
is stationary in Ti. 

Energy transfer in a light-harvesting complex 

As an example of the effect of strong interactions with a 
bath we use a model from bio-physics; energy transport 
in the Fenna-Matthews-Olson (FMO) pigment-protein 
complex, where there is thought to be significant system- 
bath entanglement and coherence [8]. As mentioned ear- 
lier, this example enables to discuss the issue of whether 
classical-correlations between system and bath can cause 
a violation of our first witness Wq (the second witness is 
not valid in this regime). 

In the methods section we impose a classical condi- 
tion based on an assumption of a class of classical states. 
States which violate this assumption possess coherences 
(either in the internal system degrees of freedom, or in 
the system-bath degrees of freedom, i.e., entanglement). 
However, to prevent classical correlations between sys- 
tem and bath from causing a false positive, the propa- 
gators Q mn (t,to) in our witness Q, which we construct 
by preparing the system in one (or more) of its states, 
must also capture the classical correlations between sys- 
tem and reservoir present at time to- In the other ex- 
amples we discuss in this work, this is trivial since the 
system and bath are always in a product state. However, 
in systems like the FMO complex we discuss here, this 
is not the case. Thus to account for these correlations 
when constructing Q mn (t, to) in a general case we must 
do the following: prepare the system-bath product state 
at t = 0, evolve to time to, and perform a measurement 
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FIG. 5: Detecting vacuum Rabi oscillations in a lossy 
cavity. Here we use the experimental data reported in Ref. 
[47] to illustrate coherence-verification using our second wit- 
ness Eq. The circular Rydberg states with principle quan- 
tum numbers 51 and 50 for transition ojq = 51.1 GHz are 
considered as the states \e) and \g), respectively. The atom- 
field coupling is ujr/2tv — 47 KHz. For a high-Q cavity with 
Q = 7 x 10 7 , the vacuum Rabi oscillation is detected by use 
of Wnmn where m = n = 1. As a comparison we also checked 
the case when the Q-factor is so low that 2lur < ujq/Q. For 
such a low-Q cavity (e.g., Q — 7 x 10 5 ), the state evolution is 
in the regime of irreversible transitions and obeys the classical 
constraint (3). Hence the value of the witness is zero. 

on the system to project it, without preserving coher- 
ence, onto one of it states n. We then evolve again, re- 
taining the post-measurement system-bath state, and de- 
duce the propagator by measuring the occupation of the 
state m at final time t. If we can do ideal projective (non- 
coherence preserving) measurements this accounts for the 
classical system-bath correlation loophole (as long as we 
can consistently prepare the t = separable system-bath 
state). If we are doing destructive or invasive measure- 
ments then we must be able to re-prepare the destroyed 
system state, at time to, on a time scale faster than the 
bath/environment dynamics. Since there is no need for 
measurements on superpositions of basis states, this pro- 
cedure can be performed without quantum tomography. 

We illustrate this with the FMO complex, a seven-site 
structure used by certain types of bacteria to transfer 
excitations from a light-harvesting antenna to a reac- 
tion center. It has been the focus of a great deal of 
attention due to experimental observation of apparent 
"quantum coherent oscillations" at both 77 K and room 
temperature. To fully capture the non-Markovian and 
non-perturbative system-bath interactions of this com- 
plex system we employ the Hierarchical equations of mo- 
tion |5j , an exact model (given a bath with a Drude 
spectral density) valid for both strong system-bath cou- 
pling and long-bath memory time. We use the param- 
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0.8^ 




FIG. 6: Detecting coherent rotations of photonic 
qubits. The first witness Wq, Eq. Q, adapted from the 
time-domain to the domain of the angles of several applied 
transformations, detects quantum coherence in almost the 
whole range of the prepared states po (see text) as a function 
of different angle settings of wave plates (4>,9): < <f> < n/4 
and < 9 < ir/2. 

eters used by Ishizaki and Fleming in Refs. [8], and 
in Figure 4 we show how this model is detected as quan- 
tum by our witness Wq , even at room temperature. We 
also show, in Figure 4c and 4d, how only partial infor- 
mation about the terms in propagator is needed to find a 
detection at small times, thus reducing the experimental 
overhead. In constructing the propagator terms for the 
sum in Eq. Q in this case we discard all coherence terms 
in the physical density matrix but retain the state of the 
bath, as in [43] . In this way we account for the state of 
the bath at time to, as discussed above. However, ac- 
counting for the classical correlations with the reservoir 
seems beyond the capability of current experiments. We 
also point out that the full witness detects coherence on 
timescales greater than to = 0.3 ps at 77 K, which is 
a much larger detection window than the Leggett-Garg 
inequality (0.035 ps) for the same parameters [44] . 

Vacuum Rabi oscillation in a lossy cavity 

We now consider a Rydberg atom placed in a single- 
mode cavity which is in resonance with an atomic tran- 
sition frequency, luq, for an adjacent pair of circular Ry- 
dberg states [IS] |e) and \g). Let us know consider the 
case when the cavity field are initially prepared in the 
excited state |e) and the vacuum state |0) p , respectively 
(denoted by |1) = |e) |0) p ). In this case, the atom-field 
state becomes |2) = \g) |l) p due to spontaneous emission 
and then periodically oscillates between the states |e) |0) p 
and \g) |l) p at the vacuum Rabi frequency u>r. If the field 
irreversibly decays due to photon loss out of the cav- 



ity, the atom-field stochastically evolves to |3) = \g) \0) p 
from 1 2). Summarizing the above, the time evolution of 
the atom-field state p can be described by the following 
master equation[46] 

-t,P= -tWjcp] - ~(afap + pa) a) + napa) (9) 
at n 2 

where Hjc = ^p-(a<r + + a^er_) is the interaction Hamil- 
tonian of the system. Here n — loq/Q, and Q is the 
quality factor of the cavity. 

We now use our second witness to detect the vacuum- 
Rabi oscillation between the atom and cavity field states. 
Here we choose the time-domain set as Tf {io[j,fc] = 
{k+j- l)t , ip,*] = (k + j - l)i<, + r\k = 1, 2, 3} for 
j = 1,2. Figure 5 shows the value of the witness for 
vacuum-Rabi oscillations in a high-Q cavity. Using the 
experimental parameters from [47] , where 2lor 3> loq/Q, 
the damped coherent oscillations of the atom-cavity state 
are detected as quantum by our second witness, shown 
in Fig. 5a. In comparison, for a low-Q cavity, where 
2ur < uio/Q, irreversible spontaneous emission out of 
the cavity will dominate the state evolution. The value 
of the witness Wfi mfl is zero for this case. The measure- 
ments on atom states we require to construct the witness 
are experimentally available by using ficld-ionization de- 
tectors [45 for selecting atom states |e) and \g). 

Coherent rotations of photonic quantum bits 

Photon polarization states \H) (horizontal) and \V) 
(vertical) have been widely used to achieve linear optical 
quantum information processing, quantum communica- 
tion, and quantum metrology [3][1H1[1S]- As a qubit, po- 
larization states can be coherently manipulated by half- 
wave plates (HWP) and quarter- wave plates (QWP). Ar- 
bitrary qubit rotations can be performed by using these 
linear optics elements. Here we will use our first quantum 
witness Wq to detect the quantum coherence of polar- 
ization states created by these rotations. The transfor- 
mations of HWP and QWP can be represented by the 
following[50 : 

H wp {<P) = cos{2<f>)(\H){H\ - \V){V\) 

- S m(m\H)(V\ + \V)(H\), (10) 

Q™ P ( ) = 4? [(*' - cos ( 20 )) \ H )( H \ + (* + cos(20)) |VX^| 
v2 

+ sm(2e)(\H)(V\ + \V)(H\)} (11) 

As a concrete example, one can set a HWP at <fi = tt/8 
to create a photonic Hadamard gate H wp (w/8). 

To detect the coherent rotations created by R(<p, 9) = 
Q wp (9)H wp ((j)), we use the first quantum witness to probe 
the coherence between states \H) and \V). While the 
witness is originally constructed in the time domain, it 



can be rephrased in terms of the settings (4>, 0) . Assum- 
ing that both the wave plates are perfect and there is no 
photon loss in the birefringent crystals of the wave plates, 
we have the following correspondences: 



(QM 0)) = tr[|m)(m| R(<j>, e)p R\4>, 0% (12) 



and 



<W<M) = \{m\R{ct>,6)\n) 



(13) 



where po is some initial state created by R. Here m — H 
and n — V denote the different measurement basis for the 
horizontal and vertical polarizations. In this example, we 
set the initial state as po = R)(4>,6) \m)(m\ R(cf>,9) and 
then the witness becomes 



1 [10 + 2 cos(40) + 2 cos(40 

16 



- cos 



COS 



(14) 



Figure 6 shows this quantum witness for different pre- 
pared states po, as a function of the angles and 4>. 

The usual approach to strictly probe the coherent su- 
perposition of states \H) and \V) is via quantum state 
tomography [50] . Compared to such tomographic mea- 
surements on single qubit states, which require three local 
measurement settings, only one setting of a local mea- 
surement is now sufficient to implement our first witness. 



Discussion 

In summary, we have formulated a set of quantum 
witnesses that allow the efficient detection of quantum 
coherence, without the restriction of non-invasive mea- 
surements. Compared to some of the existing methods, 
such as the Leggett-Garg inequality or employing gen- 
eral quantum tomography, our approach can drastically 
reduces the overhead and complexity of unambiguous ex- 
perimental detection of quantum phenomena, and has 
a larger detection window. As illustrated by the five 
physical examples, these witnesses are robust and can be 
readily used to explore the presence of quantum coher- 
ence in a wide-range of complex systems, e.g., transport 
in nanostructures, biological systems, and perhaps even 
large-arrays of qubits used in adiabatic quantum com- 
puting 51 . After this paper went to press, we became 
aware of this preprint [52 , which has related results. 
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Methods 

Proof of equation (1). The quantum two-time state- 
state correlation (Q m {t)Qn(to)} q is defined by [51] : 

(Qm(t)Qn(to))Q 
= t*SR[pSR(0)Qm(t)Qn(to)] 

= tv s {Q m (0)^ R [U(T)p SR (t )Q n (0)U^T)}}, (15) 

where psn(t ) is the system-reservoir state and U(t) is 
the system- reservoir evolution operator for t = t — t$. If 
psR.{to) is a classical state with no coherent components, 
then we have 



PSR(*o)Qn(0) = Pn(t O )Qn(0)R{t O ) 



(16) 



where p n (to) is the probability of measuring the system 
state n at time to for the classical mixture psn(to), and 
R(to) is the reservoir state at time to (which in principle 
depends on the measurement result Q n if the system and 
reservoir are classically correlated, i.e., are separable but 
in a mixture of product states). Then we have 

(Qm(t)Qn(to)} Q 
= Pn (*o)tr s { Q m (0)tifl [U{r)Q n (0)R(t o )tf (r)] } . 

The term describing the system's evolution 
tT R [U(T)Q n {0)R(t )U\T)] can be described by the 
operator-sum representation [2Tl [30] : 



tT R [U(r)Q n (0)R(t)UHr)} = £ ^(r)Q n (0)£;t(r), 

i 

where Ej(r) — ^ fe ^/p r k (rj \ U(t) jr^). The the reservoir 
state is assumed to be R(t ) = J2kPrk\rk)(rk\- Hence 
the correlation (Q m (t)Qn(to)} q for the system-reservoir 
classical mixture at the time t n is 

(Qm(t)Qn(t )) Q 

= p n (to)tr s {Q m (0)^^(r)Q n (0)E](r)} 



Pn(*o)I>^(Mo) 
3 

Pn{to)^mn{t, *0 ) : 



(17) 
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where Q mn (Mo) : = J2j^mn(t,t ) is the propagator, 
i.e., the probability of finding the state m at the time t 
when the state at an earlier time t$ is initialized at n. 

The Hierarchy model for FMO. The Hierarchy 
model was originally developed by Tanimura and Kubo 
[41"] , and has been applied extensively to light-harvesting 
complexes [8]. We will not give a full description 
here, but will just summarize the main equation and 
parameters. It is always assumed that at t = the 
system and bath are separable p(0) — ps{0) <8> Ps(0), 
and that the bath is in a thermal equilibrium state 



MO) 



'/Tr 



(3 = 1/K B T. The bath 



is assumed to have a Drude spectral density 



2 ' 



(18) 



where 7j - is the "Drude decay constant" and each site j is 
assumed to have its own independent bath. In addition, 
Xj is the reorganisation energy, and is proportional to the 
system-bath coupling strength. The correlation function 
for the bath is then given by, 



Cj — S Cj,m G X P i~^j,mt) 



(19) 



m=0 



where /i^o = Jj, & n d Hj,m = 27rm/fi./3 when m > 1. The 
coefficients are 



c i,o = 7jAj (cot(/3ft7j/2) - i) /h 



and 



4A 



Cj',m>l 



2 



7? 



(20) 



(21) 



Under these assumptions, the Hierarchy equations of 
motion are given by, 



N K 
j — l m—0 



N K 



>Pn+ 



j=1 rn— 

- *X] XI n l' m ( CmQ l P n- m - C mPn- m Qi) ■ ( 22 ) 
j — l m—0 

The operator Qj = is the projector on the 

site j, and for FMO there are seven sites, thus 
N = 7. The Liouvillian L describes the Hamil- 
tonian evolution of the FMO complex. The label 
n is a set of non-negative integers uniquely spec- 
ifying each equation; n = {n\, n^, TI3, rijv} = 
{{n w ,nii,..,niK},--,{nNo,nNi,--,n NK }}- The density 
matrix labelled by n = = {{0,0,0....}} refers to the 
system density matrix, and all others are non-physical 
density matrices, termed "auxiliary density matrices" . 
The density matrices in the equation labelled by n 



indicate that that density matrix is the one defined by 
increasing or decreasing the integer in the label n, at the 
position defined by j and m, by 1. 

The hierarchy equations must be truncated, which is 
typically done by truncating the largest total number of 
terms in a label N c = m n.^ m . This value is termed 
the tier of the hierarchy. The choice of N c should be de- 
termined by checking the convergence of the system dy- 
namics. Here we also use the "Ishizaki- Tanimura bound- 
ary condition" 42 ; 

N 00 

iiT- BC = -E E —[Qi>[Qi,P«\]' (23) 
This can be summed analytical, which for K — gives, 



rn — 1 J: 



4Aj 



{l- 7j fi[cot(7 J -^8/2)] J 8/2}. (24) 
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